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Abstract 

We study the ground state energy for a system of two hydrogen atoms coupled to 
the quantized Maxwell field in the limit a — > together with the relative distance 
between the atoms increasing as a~ 7 i?, 7 > 0. In particular we determine explicitly 
the crossover function from the R~ e van der Waals potential to the R~ 7 retarded van 
der Waals potential, which takes place at scale a~ 2 R. 

1 Introduction 

In a now very famous contribution, Casimir and Polder [I] investigate the ground 
state energy, E(R), of a system of two hydrogen atoms for which the two immobile 
nuclei are separated by a distance R and the two spinless electrons are coupled 
through the quantized Maxwell field according to non-relativistic QED. In the 
approximation where the quantum fluctuations of the Maxwell field are ignored, 
only the electrostatic Coulomb interaction remains. In this case E(R) — E(oo) ~ 
— i?~ 6 , the ubiquitous van der Waals potential, which has been discovered on 
thermodynamic grounds way before the advent of quantum mechanics. The i?~ 6 
behavior is well understood quantum mechanically [2] and has been proved in 
great generality by Lieb and Thirring [3j. Casimir and Polder use fourth order 
perturbation theory to argue that because of retardation effects the true asymp- 
totic behavior is in fact E(R) — E(oo) « —R~ 7 for large R. Their argument has 
been reworked many times and extended to arbitrary atoms and molecules, see 
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for example [U El El El E] • It is generally agreed that for two neutral molecules 
A, B it holds 

23 

E(R) - E(oo) = -—a A a B R~ 7 (1.1) 

4:71 

for large R. Here a a, &b, are the electric dipole moments of molecule A, B. The 
numerical prefactor is universal (23/47T is the value in Gaussian units). 

(II. ip is based on perturbation theory and thus holds only for small coupling. 
With improved experimental techniques, there has been a renewed interest to ex- 
plore a wider regime. One still finds the R~ 7 power law, but the prefactor is now 
a bilinear form in the electric and magnetic dipole moments. To be consistent, 
in principle, these moments have to be computed for the single molecule in isola- 
tion but still coupled to its own quantized radiation field. All atomic/molecular 
properties appear through the electric and magnetic dipole moments. As in (11.11) . 
the remaining coefficients are universal. In particular the coefficient 23/47T for 
electric dipole-electric dipole contribution persists. To mention only the most 
recent work: in [TT] IT2] the retarded van der Waals potential is computed in the 
framework of macroscopic QED. The approach in [13] is based on the standard 
non-relativistic QED hamiltonian, but uses the representation in terms of a func- 
tional integral. Conceptually this has the advantage that E(oo) is subtracted 
without error and that 1/R turns into a small parameter explicitly showing up 
in the action. Thus 1/R can be used as an expansion parameter, which is more 
physical than the conventional coupling strength to the Maxwell field. 

In the framework of non-relativistic QED the existence of a ground state, for 
arbitrary R and coupling strength, has been established in the break through 
contribution of Griesemer, Lieb, and Loss [H]. To determine the leading, large 
R asymptotics of E(R) seems to be a difficult problem, even for small, but fixed, 
coupling. In view of this situation we develop here a novel approach somewhat 
closer in spirit to the original Casimir-Polder considerations. As interaction 
strength we use the fine structure constant a and regard the ground state en- 
ergy, E(R) = E a (R), as depending both on R and a. We then study the limit 
of small coupling with an approximately adjusted scale of R, more precisely we 
consider the limit 

E a (a^R)-E a {oo), 7 >0, (1.2) 

as a — > 0. Depending on the value 7 distinct features of E a (R) will become 
visible. In particular, we will find an explicit formula for the crossover from i?~ 6 
to R~ 7 , which occurs at scale a~ 2 . 

In Section 2 we define the hamiltonians and provide an overview on the de- 
pendence on 7. There are two special values. At 7 = 1 one crosses from core 
dominated behavior to the —R~ 6 van der Waals and at 7 = 2 one crosses from 
—R~ 6 to —R~ 7 . The corresponding crossover function is computed explicitly and 
seems to be novel. Sections 3 and 4 provide proofs and point out open problems. 
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2 Hamiltonians and main results 

Let us first consider a single hydrogen atom with an infinitely heavy nucleus 
located at the origin. The nucleus has charge e, e > 0, the electron has charge 
— e. We will use units in which h = 1, c = 1, and the bare mass of the electron 
m = 1. Then the fine-structure constant is a = e 2 /4ir. Let x,p be position and 
momentum of the spinless electron. Then the non-relativistic QED hamiltonian 
for this system reads 

H 1>a = \:{p- eA(x)) 2 : -e 2 V v (x) + H { . (2.1) 

The electrons and the nuclei are assumed to have the same prescribed charge 
distribution tp with the following properties: <p is normalized, J dxcp(x) = 1, 
rotation invariant, <p(x) = <p Ta d(\ x \)i an d of rapid decrease. Denoting Fourier 
transform by <p, the potential V v is the smeared Coulomb potential 

V v (x)= [ dk\fi(k)\ 2 \k\- 2 e-' lk - x . (2.2) 

A(x) is the quantized vector potential and iff is the field energy. These are defined 
through a two-component Bose field a(k, A), k G M 3 , A = 1, 2, with commutation 
relation 

[a(k, A), a(k', A')*] = S xy 5(k - k') . (2.3) 

Explicitly 



Hf=y2 I dku(k)a{k,X)*a{k,X) 



(2.4) 



with dispersion relation 



and 



u{k) = \k\ (2.5) 



A(x)=J2 [ dkfi(k)^=^=e(k,\)(e ikx a(k,\)+e-' lk - x a(k,\y) 

= A + (x) + A~(x) (2.6) 

with the standard dreibein e(k, l),e(k, 2), k = k/\k\. : • : denotes normal ordering, 
which will be of use later on. Thus the Hilbert space for H is 



n = L 2 (Rl)®$, (2.7) 
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where $ is the bosonic Fock space over L 2 (M 3 ) (g> C 2 . From the quantization of 
the classical system of charges coupled to the Maxwell field it follows that for the 
smearing of A(x) and of Vp the same charge distribution has to be used. We refer 
to [15] for details. The ground state energy of Hi a is denoted by Ei jCc . 

To investigate the van der Waals potential one considers two hydrogen atoms, 
one located at and the other at r = (0,0, R),R > 0. It will be convenient to 
define the position of the second electron relative to r. Then X\, X2+r are positions 
and Pi,P2 the momenta of the two electrons. The two-electron hamiltonian reads 

Hr =\ : (pi - eA( Xl )) 2 : -e 2 V v {x x ) + \ :(p 2 - eA(x 2 + r)) 2 ; -e 2 V v (x 2 ) 

+ H { + e 2 V R (x 1 ,x 2 ) (2.8) 

with the interaction potential 

V R (xi,x 2 ) = -V ip (x 1 - r) - V ip (x 2 + r) + V^(r) + V v (r + x 2 - x x ) 

dk |£(A;)| 2 e ifc ' r |fc|- 2 (l - e - ifc ' xi )(l - e ik ' X2 ) . (2.9) 



Hr acts on the Hilbert space L 2 (IR^ i ) ® L 2 (M^. 2 ) <g> H R has a unique ground 
state with energy E a (R). It is known that lim^oc E a (R) = 1E\^ a . 

We plan to study E^ioT^R) in the limit of small a and consider first the 
hydrogen atom. In the limit a — > the Bohr radius is order a -1 and the energy is 
order —a 2 . Hence it is convenient to switch to atomic coordinates which amounts 
to the unitary transformation 

U*a(k,X)U = a- 3 a(a- 2 k,X) , U*xU = a~ x x , U*pU = ap, 

U*XjU = a~ 1 x j , U*p j U = ap j , j = l,2. (2.10) 

Then 



with 



U*H 1>a U = a 2 (± :(p- V4^a 3/2 A a (x)) 2 : -V a (x) + H { ) (2.11) 



V a (x)=4ir dk\if(ak)\ 2 \k\- 2 e- lk - x (2.12) 



and 



A a (x) = V / dk0(a 2 k)—^e(k,\)(e iakx a(k,\) + e~ iak - x a(k, X)*) 

A tiV R3 V 2 \ k \ 

= A + a {x) + A-{x) . (2.13) 
We note that 

a 3 {A+(x),A-(x)] = a 3 [ dk ^(a 2 ^)! 2 ^!" 1 = ^(a" 1 ) . (2.14) 
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Thus normal ordering is introduced to subtract these more singular contributions. 

Correspondingly the atomic scale hamiltonian for two hydrogen atoms sepa- 
rated by a distance a~ x R reads 



U*H a -, R U =a 2 (± :( Pl - yfazo?' 2 A a {x x )) 2 : + \:{p 2 - V^a 3 / 2 A a (x 2 + r)) 2 : 
- V a { Xl ) - V a {x 2 ) + V a>R (x 1 ,x 2 ) + H { ) . (2.15) 



2.1 The scale < 7 < 1 

For 7 = 1, instead of considering merely the ground state energy, a more complete 
picture would be the strong convergence of resolvents. For the smeared Coulomb 
potentials it holds 

lim sup I V a (x) — I re I 1 1 = , (2.16) 

Q— >o x 

lim sup \ V a>R (xi,x 2 ) 

-(-\x x + r\~ 1 - \x 2 + r| _1 + iT 1 + \r + x 2 = 0. (2.17) 

Thus the issue of strong resolvent convergence is reduced to the study of the free 
particle hamiltonian 

Ti, a = \ :{p - V^a 3/2 A a (x)) 2 : +H { (2.18) 

and correspondingly for two free electrons, with hamiltonian denoted by T 2y0l . 
Note that the norm of the coupling function in (I2.18P diverges as a" 1 / 2 . Thus 
the limit a — > is singular. On the other hand the recent estimate [16] of the 
ground state energy E® a of T^ a establishes that E® a = —a + a^a + 0(a 2 ). For 
us only the coefficient ao is of interest, which is given by 

a = (2tt) 2 <A+(0) • Af(0)n, (\Pf + H.y'AfiO) ■ A+{0)n) . (2.19) 

Here Q denotes the Fock vacuum and Pf is the field momentum, 

Pf=J2 [ dkka(k,\)*a(k,X). (2.20) 
With this information one arrives at 

Conjecture 2.1 In the sense of strong convergence of resolvents, 

lim T 1>Q = \p 2 + Ht - a , (2.21) 

a— >-0 

lim T 2 , a = \ V \ + \ V \ + H f - 2a . (2.22) 
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To our surprise, this limit has apparently never been investigated. In Section 
4 we provide some arguments towards the validity of Conjecture 12.11 If it holds, 
then by (12 .16]) and (12.17P we conclude that 

lim a~ 2 U*H x JJ = -p 2 - ■-- + H { - a 

a^O 2 \x\ 

= H hy + H i -a , (2.23) 

and 

1 1 

v -2tt*zj tt 2 i 1 2 I i-l I i-l I i-l I i i-l 

lima U H a ~i R U = -p l + -p 2 - \xi\ - \x 2 \ —\X\ — r\ - \x 2 + r\ 

a— i-0 2 2 

+FT 1 + \r + x 2 - xil" 1 + H { - 2a 
= H 2! r + H { - 2a (2.24) 

in the sense of strong resolvent convergence. Denoting the ground state energy 
of H 2t R by E 2 r, in particular it holds 

7 = 1: lim a- 2 E a (a~ l R) = E 2<R - 2a . (2.25) 

a— >0 

Note that 

\im(E 2tR -R- 1 ) = E hc (2.26) 

R^O 

with Ey ic the ground state of the helium atom, while 

lim R\E 2A - 2E hy ) = -a V w , (2.27) 

R— s>oo 

where ayw is the strength of the van der Waals potential, 

poo 

a vw = 6/ IK^o, x-e-^-^xV'o)! 2 , (2.28) 
Jo 

with H^ytpQ = Ehyip Q , Eh y = —1/2. Thus we conclude that on the distance scale 
a^R the energy a 2 E 2jR describes the crossover to the —R~ 6 potential. 
For completeness we list the even smaller distance scales, 

7 = 0: E a (R) = aV a (R) + a 2 (E hc -2a ) , (2.29) 

0<7<1: E a (a~ J R) ^ a 1+1 R' 1 + a 2 (E hc - 2a ) . (2.30) 



2.2 The scale 7 > 1 

To go beyond the distance scale a~ l R is a more difficult problem and we have 
only partial results. We expect that the range 1 < 7 < 2 is dominated by the 
van der Waals potential, i.e. 

1 < 7 < 2 : E a (a~ J R) a^~ 4 a vw R- 6 + 2E ha . (2.31) 
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The retardation of the van der Waals potential first appears at scale a 2 . 
More precisely 

7 = 2: E a (a~ 2 R) = -a 8 h co (R) + 2E 1>a , (2.32) 

where the crossover function h co is defined by 
h co (R) 

POO 

= 7T- 1 / du(l(^x ■ (H hy - E hy )((H hy - E hy ) 2 + (u/2) 2 y 1 x^))\- Ru 
Jo 

x |2- 3 i?- 2 u 4 + 2" 1 ir 3 u 3 + 5 • 2 _1 iTV + 6iT 5 u + 6/T 6 } . (2.33) 
At small distances 

h co (R) 2* a vw R~ 6 , as R -> , (2.34) 

and at large distances 

h co (R) ^ a CP R~ 7 , as R ->■ oo . (2.35) 

Here the strength of the retarded van der Waals potential is 

23 2 9 

a C P = ^(«hy) 2 , «hy = - (iphy, X • (-FZhy ~ -^hy)" 1 ^^) = ~ • (2.36) 

We conclude that at scale a~ 2 R the ground state energy crosses from the van der 
Waals potential to the retarded one as specified by h co . 

At even larger scales one expects the exact power low i?~ 7 , 

7 > 2 : E a (a~^R) ^ -a^acpR' 7 + 2E 1>a . (2.37) 

The hydrogen atom ground state has been estimated up to 0(a 5 log a -1 ) [T7j 
based on a method originally devised by Hainzl and Seiringer [18]. It is rather 
natural to use similar methods for the case of two hydrogen atoms. If more 
modestly we strive for a precision of order a 3 , then the scale will be limited to 
a~ 6 / 5 , unless there is a more direct way to accomplish the subtraction. At scale 
a~ 2 the first term is order a 8 . There is no hope to control E\ M with such a 
precision and one has to look for alternative schemes. 



3 The 7 = 2 crossover function 

Our main goal is to derive the crossover function of (12.331) . The starting point 
is the functional integral representation of E a (R) — E a (oo), see [13]. In this 
paper we consider only the second cumulant of the action, assuming that higher 
cumulants decay at least as -R -8 for R — > oo. The functional integration is defined 
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with respect to two independent ground state processes for the hamiltonian H of 
(12. ip . On the basis of our conjecture, for small coupling we replace H by 

a 2 (H a + H { -a ), H a = \p 2 - V a (x) . (3.1) 

Denoting both approximations by [• • • ] cu we have 

[E a (R) - E a (oo)} CVL = -2(2ira) 2 (l 2 (R) + I 3 (R) + I A (R)) , (3.2) 

where the coefficients are given in (39) resp. (47), (52) of [TB] with the under- 
standing that H is replaced by a 2 {H a + iff — do) • (I3.2p should be regarded as the 
definition of the left hand side. Required is an asymptotic analysis of the integrals 
I 2 , I3, and I4. In fact, in the scalings of 02.31]) . (I2.32|) . (12. 3 7p . the contribution 
of I 2 and ^3 vanish as a — > 0. The details are lengthy and will not be recorded 
here. In the following we focus only on the relevant contribution 7 4 . 

For notational simplicity, we replace H a — inf spec(H a ) by H a in the remainder 
of this section. Then the ground state ip a is defined through H a ip a = 0. Similarly 
Hhy — Ehy is replaced by H\ iy . Hence ifhyf/v = 0. According to (52) of [13], I±{R) 
is defined by 



h(R) = a 6 [ dtidt 2 dt 3 y"/ d^dA; 2 |^(a 2 A;i)| 2 |^(a 2 fc 2 )| 2 e i(fcl 

JR :i s \ JR 6 



+k2)-ra 2 



xiwiw 2 e- Wl|<1+ * 2+ * 3| e-^ 21 * 31 

x(^, ( £l • x)H a e- ik ^ xa H~ 2 e-^ e- ik * xa H a (e 2 ■ x)^ a ) 
x <^ Q , (e 1 ■ x)Hj k ^ xa H- 2 e~^ H " e ik2Xa H a (e 2 ■ x)^ a ) , (3.3) 
where now the inner product is in L 2 (IR^). 

Proposition 3.1 Assume that the smearing function ip is radial, continuous, 
and of compact support. 

(i) Let I < 7 < 2. Then 

lim a 6 - 67 (47r)- 2 / 4 (a- 7 i?) = a YW R" 6 . (3.4) 

a-¥0 



(ii) Let 7 = 2. Then 
(Hi) Let 7 > 2. Then 



lim a- 6 h(a~ 2 R) = h co (R) . (3.5) 

Q— 5-0 



23 

lim a 8 - 77 (4vr)- 2 J 4 («- 7 J R) = -^{a hy ) 2 R- 7 . (3.6) 

a— s-0 4:71 
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Using f)3.2p the proposition supports the claims of Section 2. 

Proof: For better readability we subdivide our proof into several steps. But 
before we remark that, within the current proof, compact support of ip is required. 

Step 1 (Rewriting). We scale kj ~> ay~ 2 kj and t 3 ~> a 2_7 t 3 . Then 
dtidt 2 dt 3 V [ dhdhlipia'h^lipia^h)] 2 



a 77 " 8 



xe i(fcl+fc2) - r ±u;iu; 2 e 



a 1 2 uii\t 1 +t 2 +a 2 7 t3|g-^2|i3 
" / ~ 1 ki-x tt— 2 —\ti\H a —ia~ l ~ 1 k2-x 



x(^ a ,(e 1 -x)H a e- ia1 ' k - x H~ 2 e~^ H " e" iQ7 " k2X H a {e 2 ■ x)^ a ) 
x(^, (e 1 ■ x)H a e iaJ - lk - x H~ 2 e-^ H - ^ lk ^H a {e 2 ■ x)^ a ) . (3.7) 
Let us note the following equality, 

/ dtidt 2 dt 3 e-^^'^^+^+^'^^e-^^^e-^ 1 * 1 ^^^ 1 * 21 

(0 x_i 2^ f j ^2 2Ai 2A 2 

~~ I J « y M aM Wl 2 + a 4 - 27 M 2 'a; 2 + a 4 - 27 M 2 'A 2 + M 2 'A 2 + M 2 ' 1 j 

which is proven by using the Fourier transform 



(^)- 1 /d«e-^- 5 ^- 5 = e-«W. (3.9) 
Viewing Ai and A 2 as spectral parameters for H a , one arrives at 
h(a~ J R) = a^- Q {2TxY 1 [ du [ dk 1 dk 2 \0(a^k 1 )\ 2 \0(a^k 2 )\ 2 

h 2 h 2 

A; 2 + a 4 ~ 27 u 2 A; 2 + a 4 - 27 n 2 
x2<V' Q , (ei • x)H a e- lk ^ a H~ l (H 2 a + u 2 )~ l e- ik *- xa H a (e 2 ■ x)^ a ) 

x2<^ a , ( £l • x)H a ^ xa H-\H 2 a + u 2 )' 1 e ik *- xa H a (e 2 ■ x)$ a ) . (3.10) 

Next we use 

H a (e ■ x)ip a = [H a ,e- x}ip a = is ■ pip a (3.11) 

and also introduce the integral kernel of (1 — P a )H~ l (H^ + u 2 ) -1 as 

K u>a {x,x') = (x\{t-P a )H-\H 2 a + u 2 y l \x') (3.12) 

with P a the projection onto ip a . Note that ((5 ■ p)ip a ,e lk ' x ip a ) = 0, which allows 
one to insert 11 — P a . Since H a has a spectral gap, uniformly in a, K UjCt is bounded 
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and (0, K Uj0l (j)) = C(<f),<j))u 2 for u — > oo. With this notation we arrive at the 
starting representation of I4, 



■ )\ 2 



h(a^R) = a 67 - 6 (2vr)- 1 4 f du [ dk 1 dk 2 \v(a" t k 1 )\ 2 \cp(a" t k 2 

Jr Jr 6 

xk 2 (k 2 + a^u 2 Y x k\{kl + c^V)- 1 

x ( V / dxdx'dydyX,^, x')# u , a (z/, y ')e ifcl '^ +Q7 " 1 ^- a: ))e ife2 -( r+Q7 " 1 ^'- a: '» 
x(ei -px)^! -Py)^ -2v)( £ 2 • Py')^ a {x)^ a {y)^ a {x')^ a {y') \ . (3.13) 



Step 2 (Error estimate for vanishing phase). We deal with the set on which the 
phase in (13. 13p is close to and define 

K, a = {(x,y) e M 3 x R 3 \r + a l - 1 (y-x)\>\R], (3.14) 

correspondingly A' ra with x, y replaced by x',y'. 7 4 is J 4 from (I3.13P with the 
integration restricted to A rjQ x A' ra . The error term equals If Tor — I4 — J4. 
We use the Cauchy-Schwarz inequality inside (I3.13P and perform the u,ki,k 2 
integrations. Since K ua is bounded, this yields 



r error/ - 7 



R)\<Ca~^- % [ dM^IAriHfel 2 !^!)! 2 ]^^)! 2 

x ( [ dxd2/|V^(x)| 2 |V^(y)| 2 ) 1/2 . (3.15) 



The ground state ip a has the exponential decay. Therefore 

\Il" OT {a- y R)\ < Ca-^e-^ 1 " 1 , (3.16) 

which tends to as a — > 0. In the remainder we will study I4. 

Step 3 (Angular integration). The k%, k 2 integrations are done in spherical coor- 
dinates setting dkj = w 2 dwjdflj, j = 1, 2. Let Q{k) = t— \k){k\ be the transverse 
projection. Then the angular part reads 

J dfi 1 e ifcl ' ai g(A; 1 )® j d{l 2 e ik2 - a2 Q(k 2 ) (3.17) 

with 

ai = r + oP~ l (y-x), a 2 = r + a y -\y'-x'). (3.18) 
We omit the index and compute j itt ^Q(k) for general a. 
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Let O a be an orthogonal transformation in M. 3 such that O a a = \a\e3, where 
e 3 = (0,0, 1) T . Then 

J dQe ika Q(k) = J dfie ife ° Qa O; 1 g(fc)O a 

/•7T 

= 2n d7?siiii?e i|fc||o|oOB *0- 1 5(7?)O a , (3.19) 
Jo 

where ^-(i?) = S^bji^d), i,j = 1,2,3, with 

= 62(1?) = -(1 + (cos??) 2 ), 63(1?) = 1 - (cos??) 2 . (3.20) 
Integrating over d yields 

J dtte ika Q(k) = 27rO- 1 fi(|A;||a|)O a , (3.21) 

where -By(s) = Sijbj(s) with 

b 1 (s) = b 2 (s)=g(s)-g"(s), b 3 (s) = 2(g(s) + g"(s)) , g(s) = S - 1 sms. 

(3.22) 

Thus, for j = 1,2, 

J dQj e ik ^Q{kj) = 2irO-iB{\k j \\a j \)O a . . (3.23) 
Step 4 (Radial integration). The radial integrations are of the form 

\ [ dwg{oPw)w\w 2 + c^-V)" 1 f(\a\\w\) (3.24) 



with f(s) = g(s) = s~ 1 sins or g"(s) = 2s _3 sins — 2s~ 2 coss — s _1 sins. Here 
£(|fc|) = |y3(fc)| 2 and we extended g to K by reflection at 0. We introduce a new 
function p by 

p( v ) = (2tt)- 1/2 / dwg(w)e [vw . (3.25) 



Then one has 

r dvp(v) = (2n) 1/2 g(0) = (2vr)- 5/2 . (3.26) 



Let cr(|x|) = if * <p(x). Then 

p( v ) = (2tt)- 3/2 / drm(r) (3.27) 



for v > 0. Since <p is continuous and of compact support, p G C 1 (M) and p has 
has compact support. 
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We plan to use Plancherel's theorem in (I3.24p and obtain, in the sense of 
distributions, 



ciO; \a\,u) = (2tt)- 1/2 / dw e lvw w 3 (w 2 + n 2 )" 1 ^!" 1 sin(|a|w) 

Jr 

= (4|a|)- 1 (2 7 r) 1 / 2 ( - 25' (v + |a|) - u 2 sgn(t; + |a|) e~^ v+ ^ 

+ 25' (v - \a\) + u 2 sgn(v - |a|) e - |u|l "- |a|1 ) , (3.28) 

c 2 (v; \a\,u) = (2tt)- 1/2 / d«;e im «)(«; 2 + m 2 )" 1 ^!" 3 sin(|a|u;) 

= (4|a| 3 )- 1 (27r) 1 / 2 (sgn(v + |a|) e" 1 " 11 ^ 11 - sgn{v - |a|) e" |u||,, - |a|1 ) , 

(3.29) 

c 3 (v; \a\,u) = (27T)" 1 / 2 / dW^w 2 (u> 2 + m 2 )" 1 ^!" 2 cos(|a|w) 
= (A\a\ 2 )' 1 (2n) 1 / 2 (25(v+\a\) - \u\ e~^ v+ ^ 



+ 25(v-\a\) - |M|e- |u||l, - |a|1 ) (3.30) 

with the sign function sgn(t) = 1 for t > 0, sgn(t) = —1 for t < 0. Since p G C 1 , 
Plancherel's theorem yields 



di(|a|, u, a) = / dva 7 p(a 7 t>) (ci(t>; |a|, u) — c 2 (v; |a|, u) + c 3 {v] \a\, it)), 
c/2 ( I a 1 5 w > a ) = di(\a\,u,a) , 

ds(\a\, u, a) — / dv a' 1 p{a~ 1 v)2(c2{v ; |a|, u) — o 3 {v ; |a|, u)) (3.31) 



and 

Dij(\a\, u, a) = (2ii)5ijdj(\a\, u, a) . (3.32) 

We now combine all terms. The ground state ijj a is invariant under rotations. 
It is convenient to write 

ip a (x) = ipa,rBd(\x\) , W> a (x) = ip' ajI3d (\x\)x , (3.33) 

where x = x/\x\. Then (13 . 13[) becomes 

= a 67 ~ 6 (27r)~ 1 4 I du [ dxdydx'dy' 

X (x ■ Q r+a^{y-x) D (\ r + (V ~ x)\, a 2 ' 1 a)O r+al -i {y _ x) y^j 

X ■ °r+a-r-Hy'-x') D (\ r + a ^ iv' ~ x')\, a 2 ~^U, a)O r+a ^i (yl _ xl) y'^j 

X <,rad(kl)<,rad(k1)<,rad(bl)^,rad(IZ/'l)^,«(^ x')K u>ot {y, y') . (3.34) 
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Step 5 (The limit a — > 0). c±, C2, C3 contain terms proportional to 5 and 5'. Since, 
by assumption, g has compact support and since \aj\ is bounded away from zero 
on the prescribed domain of integration, these terms vanish for a sufficiently 
small. Thus only the regular terms, containing the exponential function, have 
still to be considered. 

We have to discuss the cases 1 < 7 < 2 and 7 > 2 separately. 

1 < 7 < 2. As before we use the uniform bound from K Uj0l . Therefore the terms 
proportional to u°, u, u 2 have a uniformly integrable bound in u. By dominated 
convergence only the term proportional to u° does not vanish as a — > 0. The 
term proportional to u 3 , u A are bounded as 

(1 + u 2 )' 2 {{ a 2 -'uf + (a 2 ">|) 3 ) e- Ka2 ~ 7 H < Ca 3 ^ e^ 2 " 7 ^! (3.35) 

with k > kq > uniformly in a. Thus the integral over u vanishes as a — > 0. 

We are left with the products of the u° terms. As a — > 0, the matrix 
Or+c^-i^-x) tends to the unit matrix. Thus we conclude 

lim a 6 ~^h(a~^R) 



2 



(27r)- 1 4((27r) 3 / 2 (2 J R 3 )- 1 | / dvp(v) 



x 



X ^hy,rad(kl)^y,rad(k1)< y ,rad(bl)<y,rad(b1)^n,o(a;, x')K u>0 (y, y') . (3.36) 

Using the rotational invariance of K Uj0 one arrives at 

lima 6 - 67 /^- 7 ^) = (27r)- 3 2- 1 3i?^ 6 f du(§(^ hy , x ■ H hy (Hl y + m 2 )" 1 ^)) 2 

Jr 

= (2n)- 2 2- 2 a YW R- 6 . (3.37) 

2 < 7. We substitute u by a 7 ~ 2 w. The uniform bound now results from the 
exponential terms exp[— \u\v ± \a\], using that 

/ dw|p(w)|e- |u||a7l,±|a|1 < Ce- K|u| (3.38) 
Jr 

uniformly in a, provided a is sufficiently small, since p(v) has compact support 
by the assumption. In the limit a — > one obtains a formula which has the same 
structure as in (13.361) . Only the coefficients in front of Xj,{jj and x'j,y'j are now 
different. 

For 7 = 2, one obtains 

2 



lim a- 6 h(a- 2 R) = (2tt) 2 ( [ dvp(v) 

v Jr 

x / d M (/3 2 + 3/3 2 + 3/3 3 2 -2/3 1 /3 2 -6/3 2 /33 + 2/3 1 /3 3 )2e- 2if l u l 
Jr ^ ' 

x (±{iky,x ■ (1 - P hy )H hy (Hl y + n 2 )- 1 ^^)) 2 (3.39) 
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with Pi = —R 1 u 2 , 02 = R 3 ; 03 = —R 2 \u\, which yields 

lim a- & h{a~ 2 R) = (2ny 2 2~ 1 h co (R) . (3.40) 

a— >0 

For 7 > 2, one has the same expression except that from the rescaling of du 
one picks up the factor a 7-2 and that the factor (11 — Ph y )(-Ph y + w 2 ) -1 now reads 
(11 — PhyX-^hy + (a; 7 ~ 2 it) 2 ) _1 , which is still uniformly bounded in u. Hence 

lim a 8 - 7 ~<h(a-^R) = (2tv)~ 2 2" 1 a cv R~ 7 . (3.41) 
This concludes the proof of Proposition 13.11 □ 



4 Strong resolvent convergence 

We discuss the limit a — > for a free electron coupled to the radiation field 
on the scale set by the hydrogen atom. Then the energies are of order a 2 and 
hamiltonian on that scale reads 

T 1>a = | :(p- V^a 3 ' 2 A a {x)f: +H { . (4.1) 

The coupling function in (14. ip is 

g a (k, A) = V4^a 3/2 ^(a 2 k)^=e iak - x e(k, A). (4.2) 

Note that \\g a \\ = a" 1 / 2 , which makes the limit a — > singular. 

The ultraviolet cutoff as a~ 2 corresponds to a charge distribution localized on 
the relativistic scale. If the charge distribution would have a width of the order 
of the Bohr radius, then 0(a 2 k) would have to be replaced by ip(ak). Thus it 
is natural to introduce the parameter 5 with < 5 < 1 and to define by 
(14. ip with <p(a 2 k) substituted through (p(a 2 ~ 5 k). If < 5 < 1, our arguments in 
Sections 2 and 3 would not be altered, except for ao = 0. But now the resolvent 
convergence can be established. 

Proposition 4.1 Let < 5 < 2. Then, in the sense of strong convergence of 
resolvents, 

lim T[ 5) a = \p 2 + Hi . (4.3) 

Proof: Let 

TS = §:(p- V^o?l 2 A a j{x)) 2 : +H { = T + B a , s , (4.4) 

where 

T = ±p 2 + H f , 

B a , s = -V^a 3/2 p ■ (A+fix) + A- t5 (x)) 

+2na 3 {Al s (x) ■ A^(x) + 2A^(x) ■ A",(x) + A" 5 (x) • A~ a (x)) . (4.5) 
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The coupling function of A aj s{x) is given by 

g a ,s{k, A) = <p(a 2 - 5 k)-^e iak x e(k, A) . (4.6) 

If it can be shown that 

1(0,5^0)1 <C(a) || (T + i) 1/2 0||, C(a)40asa4fl, (4.7) 

for all <p G dom(T 1 ^ 2 ), then one concludes — > \p 2 + H{ as a — > in the norm 
resolvent sense by the general theorem [TH1 Theorem VIII. 25], as based on the 
famous Nelson's argument [2D]. To prove ( 14. 7p . we apply the standard bounds 

||a(/)^|| < l^-^/IIIK^f + 1) V V|| , (4.8) 
||a(/)>|| 2 <(||/|| 2 +||^ 1 / 2 /ir)ll(^f + ]l) 1/ ^l| 2 - (4.9) 
As to A~ s (x), the bound (14. 8[) translates to 

o?' 2 \\A^{x)n ^^ll^-^IHKTo + ll) 1 / 2 ^!! (4.10) 
< O(a (1+5)/2 )||(T + l) 1/2 tP\\ . (4.11) 

Similarly A~ s ■ A~ s can be estimated as 
a^(Hf+ti)-V 2 A- >s (x)-A- 5 (x)i,\\ 

< « 3 i<^f + n)- 1/2 ii wka^ + *r 1/2 \\ m + ty 2 n 

< a\\\g a A\ 2 + \\uj- l ' 2 g a A\ 2 ) l/2 \\^l 2 g a A ||(Jf f + H) 1/2 ^ll 

< C(a 35/2 )||(/7 f + ]l) 1/2 7/;||. (4.12) 

Thus we arrive at 

a^K^p ■ A- 6 <f>)\ < O{a^ 2 )\\{T + l)W<f>f (4.13) 

and 

^A~^x) ■ A-/x)cf>)\ < a 3 \\(H { + H)^|| \\(H { + l)' 1 ^ • A'^W 

<C(a 35/2 )||(T o + ll) 1/2 0|| 2 (4.14) 

1 /2 _____ 

for each <p 6 dom(T ). Hence (14. 7p is satisfied and the assertion follows. □ 

The case 5 = is physically distinguished, but Nelson's argument fails and no 
other functional analytic method seems to be available. We devise an alternative 
approach based on functional integrals, which clearly displays that 5 = is on 
the borderline. In our context functional integration is explained in [15] , Chapter 
14, and at greater depth in [21]. The propagator e~ rTl ' a can be written as an 
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integral with respect to Brownian motion for the particle and a Gaussian space- 
time measure for the Maxwell field. It is convenient to pick for ijj the particular 
form 

V> = 0®W(/)fi. (4.15) 
Here G L 2 (M 3 ), Q is the Fock vacuum, and W(f) is the Weyl operator 

W(f) = e(°^*-^ , / G L 2 (R 3 ) ® C 2 . (4.16) 



Note that the linear span of these vectors is dense in L (Mr) <g)#. Integrating over 
the Maxwell field one arrives at the following expression 

(V,e-^> = E w (cf>{q )*<f>(q T ) e"*) , (4.17) 

where t > qt is a path in Mr~ and Ew denotes average over the Wiener measure. 
The action results from the Gaussian integration over the Maxwell field and 
consists of a sum of three pieces, 

si = stfx + + s/ 3 . (4.18) 

srf\ is the piece corresponding to / = 0, 

M = 47ra 3 [ [ dq t ■ W a {q t - q s ,t - s) dq s (4.19) 
Jo Jo 

with the photon propagator 

W a (x,t)= [ dk\0(a 2 k)\ 2 7 ^e ik - xa e~^Q(k) (4.20) 

J R 3 AbJ[ti) 

and Q(k) = i—\k){k\, the transverse projection. (14 . 1 9 j) is an iterated Ito integral. 
It avoids the diagonal {s = t} in accordance with the Wick ordering : :. reflects 
the term coming from W(f), 

s^= f dk-(l + e-^)f*(k)-Q(k)f(k). (4.21) 
Jrs lo 

Note that does not depends on q t and a. Finally the cross term reads 

^ 2 = -iV^a 3/2 [ [ dkif(a 2 k)-^e ikqta {e- u)t + e""^) dq t ■ Q(k)f(k) . 
Jo Jr3 V 2w 

(4.22) 

The cross term is small, since for the expectation Eo with respect to standard 
Brownian motion starting at go — it holds 

E (|^ 2 | 2 ) = 4vra 3 f I dhdh^k^ofk^k^k^ 1 ' 2 
Jo Jk. 6 

xe -(a a |(fci+* a ) a +|fc 1 |+|te|)t^ A . i j . Q(k 1 )Q(k 2 )f{k 2 ) 
< aTrll^u;- 1 !! 2 !!/!! 2 . (4.23) 
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Ignoring the cross term one arrives at 
(0 ® W(f)n, e-^cf) <g> W(f)Q) = E w U(q )*<P(qr) e~^) + G(a) . (4.24) 



It is now convenient to rewrite (I4.24p using that stf\ depends only on the incre- 
ments. Then 



<0® W(f)Q,e- TTl ^<f)® W(f)fl) = / dfc|0(A;)| 2 E o (e- £/l e ifc ^)e- £/3 + 0{a) , 

(4.25) 

To turn to srf\ we first note that Eo(^) = by Ito calculus. Secondly we use 
the Brownian motion scaling q t = aq t / a 2 to rewrite srf\ as 

M = 47ra / dq t ■ W^qt - q.,t - s) dq s . (4.26) 

Jo<s<t<r/a 2 

By definition W\ does not depend on a. W\ decays as it — s)~ 2 , which should 
provide enough independence for a central limit theorem hold. Thus the key 
input is that q T ) jointly converge to a Gaussian as a — > 0. One checks that 

Eo(M?r) = 0. (4.27) 

Hence the assumption is that 

lim M = £ G (4.28) 

ct-¥0 

with a centered Gaussian random variable independent of q T . To complete our 
argument we compute the variance of sd\, 

E «) 

= (47ra) 2 / / E ((dg tl • Wx{q H - q 8l ,h - Si)dq 8l ) 

^0<si<ti<r/a 2 Jo<S2<t 2 <T/a 2 

x (dq t2 ■ Wi(q t2 - q S2 ,t 2 - s 2 ) dq S2 )) 
= (47ra) 2 f E (Tr (W^qt - q„t - s) 2 )) dsdt . (4.29) 

J0<s<t<T/a 2 



Hence 



poo r 

limE «) = (4tt) 2 / dt dfc 1 dfc2^(fci)ri^(fc 2 )| 2 (2|fc 1 |2|A ;2 |)- 1 

Jo JR6 

x exp { - l - ({k x + k 2 f + \k x \ + \k2\)t}Tr(Q(kx)Q(k2)) 
= 2a r. (4.30) 
Returning to (I4.24p one concludes that 



lim (0® W(f)n,e- TTl < a (f)® W(f)Q) = / dA;|0(A;)| 2 Eo(e ifc ^)E(e ?G )e- £/3 

= (0 ® e - T «P 2 / 2 )+^- a °)0 <g> . (4.31) 

If one reintroduces the parameter S from above, then the variance vanishes 
provide 5 > 0, in accordance with Proposition 4.1. 
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